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27. (Cancelled) 

Remarks 

[1] Claims 1 -27 are pending in this application. Claims 1-27 stand rejected. Claims 2-4, 
21-23 and 27 are cancelled. As a courtesy to the Examiner, the Applicants have used the 
paragraph structure used by the Examiner to respond to the Examiners requests, 
objections and rejection. Claims 12-19 are amended to place them in condition for 
allowance based on the Examiner's statement that they are non-obvious over the prior art 
Certain other claims are amended to place them in better form in case of the necessity of 
an appeal. No new limitations are included in the amendments; so if Examiner does 
believe they are currently allowable despite Applicants 5 firm belief that they are, it is 
respectfully requested that they be entered for purposes of narrowing the issues for an 
appeal. 

Priority 

[2] The Examiner is correct that the Supplemental Application Data Sheet is in error. A 
corrected Supplemental Application Data Sheet referencing the provisional application is 
submitted with this amendment as Appendix I. 

Drawings 

[3] The Examiner has objected to the drawings on several grounds. To overcome each 
the Applicants have done the following: 
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[3-1] FIG. 6 is amended to substitute "Fail Rate" for "Fall Rate" per the originally 
submitted drawings. In FIG. 8 "CPU" is substituted for "RAM" for reference character 
10 as per the originally submitted drawings. 

[3-2] In FIG. 4 the drawing is corrected to remove the first reference character 58, 
substituting *56" per the originally submitted drawing. 

The Applicants apologize that the formal drawing contained these mistakes and believe 
the Examiner's rejection is overcome. 

Specification 

[4] The Applicants thank the Examiner for withdrawing the objections to the 
specification* 

Claim Rejection - 35 USC §112, First Paragraph 

[5] The Examiner has rejected the application under 35 USC 1 12, first paragraph, 
arguing that there is a substantial lack of teaching for the claimed invention. 

[6] In the rejection states the Examiner states "how to 'combine the individual values and 
then add hot-e input based on performance sort" have not been described in the 
specification. 11 

The Applicants disagree. As stated in the previous amendment the Monte Carlo analysis, 
which Applicants are using to Combine the individual values and then add hot-e input 
based on performance sort' is well known in the art as demonstrated by readily available 
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textbook sources. To prove this point Applicants have attached in Appendix II three 
common text book references discussing Monte Carlo analysis in engineering 
applications. For the first two references only the title pages are provided. The third 
source, Probabilistic Systems Analysis, presents a practical discussion and application of 
the Monte Carlo technique. (Note that the author subtitles it as an "Introduction to 
Probabilistic Models/ 3 ) Applicants have attached several pages from that text The 
Examiner should particularly note the flow chart on page 157. It describes the basic steps 
behind the Monte Carlo analysis, and forms the basis for Applicant's FIG- 7. Please note 
also that Figure 7-6 on p. 173, which shows the output of a particular Monte Carlo 
analysis. (See FIG. 2.) Simply put Applicants are applying a known statistical technique 
to solve a guardband problem in a novel way by accounting for the tester environment, 
tester-to-system offset, system environment, and reliability parameters (e.g, hot-e 
degradation). 

In addition, Figure 7 and the text referencing Figure 7 describes Applicants' preferred 
Monte Carlo analysis of selecting a value from each distribution type (tester environment, 
tester-to-system offset, and system environment) and then combining (i.e. adding) those 
values with a hot-e value (based on performance sort) to determine a guardband value for 
each iteration of the analysis. Applicants* Monte Carlo function, selecting a value from 
each distribution type and combining those values with a hot-e value, is comparable to 
the generic Monte Carlo function as illustrated on p. 157 of the third source, Probabilistic 
Systems Analysis, where the Monte Carlo function is depicted as "Compute yj = g(xi 
X2o)-" process is continued for n iterations in accordance with the Monte Carlo 
approach at which time the output is the distribution of guardbands as illustrated in 
Figure 2. 

Each value contributed from the tester environment distribution model contributes a 
performance value (e-g. FMAX) as described on page 3, lines 12 on. Thus, when the hot- 
e value is added to the combined distribution values for a particular iteration of the Monte 
Carlo analysis, the hot-e value added to these combined values corresponds to the 
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performance value (FMAX) as taught at page 6, lines 8-19 and in accordance with Figure 
6. 

In conclusion, the Applicants specification is in full compliance with 35 U.S.C., first 
paragraph. 

Claim Rejections - 35 UJS.C. 112, 2nd Paragraph 

[7] The Examiner rejected claims 2-3, 12-) 3 and 15-17 under 35 USC 1 12, second 
paragraph. 

[8] Applicants have cancelled claims 2 and 3 and amended the other claims to overcome 
the Examiner's rejections. Claim 1 now includes claim 2 and has an antecedent basis 
added for "product." Claim 16 now is dependent on claim 13 which contains the 
referenced step. Claim 20 is amended so that "product" has an antecedent basis in the 
preamble. 



Claim Rejections -35 U.S.C. 102 
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[10] The Examiner rejected now pending claims 1 , 6-7, 20-24 and 27 under 35 USC 
102(b) as being anticipated by Mrttle et al., U.S. Patent 5, 634,001 issued May 27,1997 
("Mittle"). 

(10-2, 3, 4] Mittle does not teach or suggest modeling the tester environment or end- 
system environment in which the microprocessor will reside. Claim 1 has been amended 
to introduce in the claim both the system on which the product is used and tester 
offset a s variables. MitUe teaches a method of determining a hot-e guardband by taking 
into consideration circuit timing data, operating voltage and case temperature as well as 
the voltage / frequency response of a microprocessor. 

These are all internal characteristic of the of the "product." ie, microprocessor or inputs 
to the microprocessor. They are not the "system" in which the product is used as used 
and denned by applicants. Nor does Mittle speak of tester to system offeet. Nowhere in 
the section sited by the Examiner is there any mention of system to tester offset 
especially in terms of what applicants clearly refers to and claims as its invention. Mittle 
restricts its teaching to the calculation of a hot-e guardband which is only one 
consideration of the present invention. 

As stated previously and repeated here, Mittle does not teach or suggest modeling of 
"system variables." Instead, Mittle only considers a nominal operating voltage and 
temperature at which a microprocessor is expected to operate when determining a hot-e 
guardband. Therefore Mittle does not teach "creating a set of distribution models 
representative of variables'* in claim 1. 
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As stated previously and repeated here, Mittle simply determines the difference expected 
in a microprocessor's performance between beginning of life (BOL) and end of life 
(EOL). Mittle does not teach or suggest that the test environment upon which a 
microprocessor is tested should be modeled to account for variations within that system- 
Hence there is no teaching of a variable for u system to tester offset" as set forth in claim 
1 . Mittle discloses circuit-level timing analysis. This is based upon a simulation of a 
particular circuit. Mittle does not teach or suggest the use of a tester to determine actual 
circuit performance. A simulation environment on a computer does not inject variations 
into guardband determinations as does an actual tester environment An actual tester 
environment introduces many variations as a result of both electrical and mechanical 
components. These components do not exist when simulating circuit performance, and 
thus there is no need to model variations that may exist in a tester environment. So there 
is no teaching of using the "test system" as a variable. 

[10-5] The arguments with regards to claim 1 render claim 6 allowable. 

[10-6] As stated previously and as preserved here, Mittle does not teach or suggest a 
sample of at least 10. In feet, there is no sample size requirement in Mittle. The part of 
Mittle cited by the Examiner discusses a histogram of propagation delays for all the paths 
of a microprocessor. There is no sampling as set forth in claim 7, just simply a graphical 
representation of the delay associated with each path of a microprocessor. 
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[1 1-7] As stated previously and as persevered here, Mittle does use a reliability wearout 
model it does not teach or suggest the use in the context of this invention; i.e., of a "set of 
distribution models representative of variables that affect the specification:' So claim 10 
is allowable in view of how it the model is used. 

[10-8] As stated previously and as preservered here, the Examiner notes that claim 1 has 
similar limitations to those in claims 20. Therefore the arguments overcoming 
anticipation for that claim applies equally well to claim 20. 

[10-9] Claim 27 is cancelled. 

[11] The Examiner also claims 1, 8-99, 20-23 and 27 under 35 USC 102(b) as being 
anticipated by Conrad, et al, Calculating Error of Measurement on High Speed 
Microprocessor Test," Proceeding of International Test Conference, October 1994, pages 
793-801 ("Conrad"). 

[11-1] Claim 1 now contains the limitations in claims 2 and 3. Thus claim 1 is allowable 
over Conrad. 

[1 1-2,3] Claims 8-9 now rely on claim 1 which in its new form is allowable over Conrad. 

In lieu of the above arguments, it is clear that neither Mittle or Conrad teach what is 
claimed by Applicants. 
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Claim Rejections - 35 U.S.C- 103 

[13] Claims 5 and 25 are rejected under 35 USC 103(a) in view of Mittle and 
Applicants assertion. 

[13-1] As stated previously and as preserved here even if it would be obvious to use a 
Monte Carlo analysis in light of Mittle, this is irrelevant because Mittle teaching is 
limited to how to determine a guardband for reliability wearout mechanisms. Mittle fails 
to suggest the use of models (especially models relating to use of the product in a system 
and tester offset that are representative of variables in a specification* See discussions 
regarding claims 1 and 20 above. As to the Monte Carlo analysis, Mittle actually teaches 
away from the use of a Monte Carlo analysis because such an analysis would inject 
unnecessary uncertainty into the guardband calculation. The techniques taught in Mittle 
have no use for a statistical sampling technique like that taught and claimed by 
Applicants. Statistical sampling may actually cause excessive variation in the hot-e 
guardband calculation, thus rendering it ineffective. Monte Carlo analysis is suited for 
applications where a probabilistic approximation to a solution is desired. Mittle teaches 
how to determine an exact guardband based upon circuit timing analysis, operating 
voltage and temperature and channel length. A probabilistic approximation is not 
desirable in this situation because nothing needs to be approximated. 
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As stated previously and preserved here. Applicants assertions regarding the variability of 
the type of statistical analysis tools in its own claimed invention has nothing to do as to 
whether it would be obvious to include such tools in Mittle. 

Since the Examiner has failed to establish that Mittle or Conrad reads on claims 5 and 25, 
let alone the claims upon which they depend, the Examiners 103 arguments with respect 
to those claims fail. 

[14] The Examiner has rejected claims 1 1 and 26 under 35 USC 1 03 as being 
unpatentable over Mittle ii\n view of Kreyszig, 44 Advanced Engineering Mathematics " 
John Wiley & Sons, 1 988, pages 1248-1 253 ("Kreyszig"). 

[14-1] With regard to the remaining 103 rejections, the combinations cited by the . 
examiner fail to suggest; three of the novel aspects of the present invention set forth in 
claims 1 1 and 26. First, there is no suggestion of the use of "distribution models 
representative of variables that affect a specification." Furthermore, there is no 
suggestion that the distribution models should be analyzed using a statistical tool. The 
cited art only deals with individual components of a guardband (reliability wearout 
mechanisms in Mittle). There is no suggestion as to how one would determine an overall 
guardband based upon reliability wearout mechanisms as well as variations within both 
the tester and system in which the product is used. Finally, there is no suggestion to 
select a final guardband based upon a statistical analysis. 

For these reasons claims 1 1 and 26 are clearly patentable over the prior art. 
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[15] Applicants thank the Examiner for finding claims 12-19 non-obvious over the prior 
art. These claims now stand alone and are in condition for allowance. 

SUMMARY AND CONCLUSION 

In view of the foregoing, withdrawal of the rejections and the allowance of the 
current pending claims is respectfully requested. If the Examiner feels that the pending 
claims could be allowed with minor changes, the Examiner is invited to telephone the 
undersigned to discuss an Examiner's Amendment 

JResnect fullv submitted, 




BY: Richard M. Kotulak 
Attorney for the Applicants 
Registration No.: 27,712 
Telephone: 802-769-4457 
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Practitioner's Docket No. filTTmO?^. n itIK1 PATENT 

IN THE UNITED STATES PATENT AND TRADEMARK OFFICE 

In re prior application of: Mark R. Bilak 

Application No.: 09 / 520,257 Group No.: 2123 

Filed: 03/07/2000 Examiner Day, Herng-Der 

For STATISTICAL GUARD BOND METHODOLOGY 



MaB Stop Patent Application 
Commissioner for Patents 
P.O. Box 1460 
Alexandria VA 22313-1450 



SUPPLEMENTAL APPLICATION DATA SHEET 
37 OF-R. § 1.76(c) 

NOTE: "Supplemental application data sheets. Supplemental application data sheets: 

(1) May be subsequently supplied prior to payment of the issue fee either to correct or update 
information in a previously submitted application data sheet or an oath or declaration under § 1.C3 
or $ 7.67, except that inventorship changes are governed by § 1.48, correspondence changes are 
governed by § 1.33(a), and citizenship changes are governed by § 1.63 or $ 137: and 

(Z) Should identify the information that is being changed (added, deleted, or modified) and therefore 
need not contain aU the previously submitted information that has not changed 9 37 § 1.76(c). 

The fbHowing information on the Application Data Sheet is changed as indicated: 



CERTIFICATION UNDER 37 CJ.R. 6$ 1.8(a) and 1.10* 

(When using Express Ma8, the Express Mali label number Is mandatory; 
Express Mali certification is optional) 

I hereby certify that, on the date shown below, this correspondence is being: 

MAILING 

□ deposited with the United States Postal Service In an envelope addressed to Commissioner for Patents, P.O. 
Box 1450, Alexandria, VA 22313-1450 

37 C«F.R. § 1J«a) 37 O.F.R. § 1.10 # 

□ with sufficient postage as ft ret claoa man. □ as "Express Man Post Office to Addressee" 

Mailing Label No. - (mandatory) 

TRANSMISSION 
facsimile transmitted to the Patent and Trademark Office, (703) . 



Date: ^Z/^ 




{type or print name of person certifying) 



- Only the date of Wing (§ 1.6) wifl be the date used in a patent temi adjustment calculation,, although the date 
on any certificate of mailing or transmission under $ 1.8 continues to bo taken into account in determining 
timeliness. See § 1.703(1). Consider "Express Mail Post Office to Addressee 9 (§ 1.10) or facsimile transmission 
(§ 1.0(d)) for the reply to be accorded the earliest possible filing date for patent term adjustment calculations. ' 

(Supplemental Application Data Sheet [4-13]— page 1 of 6) 
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BIBLIOGRAPHIC DATA 



1. □ Applicant information is being □ added □ deleted □ modified: 

NOTE: 37 CER. $ 7. 76(b)(1); "(1) Applicant information. This information Includes the name, residence, mailing 
address, and citizenship of each applicant (§ 1.41(b)), The name of each applicant must include the 
family name, and at teas? one given name without abbreviation together with any other given name or 
initial. If the applicant fa not an inventor, thfa information also includes the applicant's authority (§§ 1.42. 
1,43, and 1.47} to apply for the patent on behalf of the frnwror." 

WARNING: Inventorship changes are governed by § 1.48, 37 CFJi, § 1.76(c)(1). 

WARNING: Citizenship changes are governed by § f.63 or § 7.S7. 37 C.F.R. 1.76(c)(1). 

First applicant: 

GIVEN NAME MIDDLE INITIAL OR NAME FAMILY (OR CAST NAME) 

Citizenship 

Residence 



Second applicant, (if any) 

GIVEN NAME ' MIDDLE INITIAL OR NAME FAMILY (OR LAST NAME) 

Citizenship ^ 

Residence ; — ; 



Third applicant, Of any) 

GIVEN NAME MIDDLE INITIAL OR NAME FAMILY (0* LAST NAME) 

Citizenship . , lm — 

Residence 



Fourth applicant, if any 

GIVEN NAME " MIDDLE tNfTlAL OR NAME FAMILY (OR LAST NAME) 

Citizenship , 

Residence — . 



(Supplemental Application Data Sheet [4-1.2}— page 2 of 6) 
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Fifth applicant Of any) 

GIVEN NAME MIDDLE MJT1AL OR NAME PAM1LY (OR IAST NAM© 

Citizenship 

Residence, , — — — 



Sixth applicant, (ft any) 

GfVBtNAME MIDOLE BOTAL OR NAME FAMILY PR LAST NAME) 

Citizenship ! _ — 

Residence - — , — — 



□ Applicant knot the Inventor and applicant's authority (§§ 1.42, 1.43 and 1.47) 
to apply for the patent on behalf of the Inventor Is as follows: 



2. Correspondence Information Is being □ added □ deleted □ modified: 

NOTE: 37 CF.R § 1J6{b)ft: m (9 Correspondence information. This Information Includes the correspondence 
address, which may be indicated by reference to a- customer number, to which correspondence Is to 
bo Greeted pee $ J,33ftp.- 

WARNING: Convspondence changes ore governed by $ 1.33(b) 

Correspondence for this application should be addressed as follows: 
Name: 

Address: 



□ Customer No.: 



(Supplemental Appecation Data Sheet (4-1-3 ) papo 3 of 6} 
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3. Application information is being □ added □ deleted □ modified: 

NOTE: 37 GF.A § 1.76(b)(3): 'Application information. This information includes the title of tfie Invention, a 
suggested classification, by class and subclass, the Technology Center to wNch the subject matter of 
the invention is assigned, the total number of drawing sheets, a suggested drawing figure for publication 
Qn a nonprovfsional application), any docket number assigned to the application, the type of application 
(e.g., utility, plant, design, reissue, provisional), whether the application discloses any significant part 
of the subject matter of an application under a secrecy order pursuant to § 5.2 of this chapter (see 
§ 5.2(c)), and, for plant applications, the Latin name of the genus and species of the plant claimed, 
as well as the variety denomination. The suggested classification and Technology Center information 
should be supplied for provisional applications whether or not claims are present, tf claims are not present 
in a provisional application, the suggested classification and Technology Center should be based upon 
the disclosure. " 

Title of Invention: 

Docket number assigned to this application: 
Suggested Classification: Class: 

Subclass: 

Technology Center to which subject matter Is assigned: 

NOTE: The suggested classification and Technology Center Information should be supplied for provisional 
applications whether or not claims are present tf claims are not present in a provisional application, 
the suggested classification and Technology Canter should be based upon the disclosure.' 37 CF.R. 
§ 1J6(b)m. 

Total number of drawing sheets: 
Type of application: 

□ utility 

□ application Is to be published 

Suggested drawing figure for publication: ^ 

□ application is not to be published 

□ plant 

□ Latin names of the genus , , 

species 

of plant being claimed* 

□ design 

□ reissue 

□ provisional 
Secrecy order under § 5.2: 

This application 

□ does not disclose 

□ discloses a significant part of the 

subject matter of an application which is under a secrecy order pursuant to § 5.2. 

(Supplemental Application Date Sheet TF5>[4-1 .2) — page 4 of 6) 



PAGE 19/25 * RCVD AT 1QM /2004 7:12:20 PM [Eastern Daylight Time] * SVR:USPT0-EFXRF-1/3 * DNIS:8729306 * CSID:8027698938 * DURATION (mm-ss):08-34 



OCT 04 2004 19=35 FR BTU IPLAW 



8027698938 TO 917038729306 P. 19/24 



4. Representative information is being □ added □ deleted □ modified: 
WOT& 37 C.F.R, § 1.76(b)(4) states; "Representative information. This information includes the registration 
number of each practitioner h&ving a power of attorney or authorization of agent in the application 
(preferably by reference to a customer number). Providing this information in the application data sheet 
does not constitute a power of attorney or authorization of agent in the application (see § 1 .34(b]). - 
(Emphasis added). 

The following have a power of attorney or authorization of agent In this application; 

Name of attorney (agent): — — — 

Address: ~ — — — 



Customer No.: _ — 

5. Domestic Priority information is being □ added □ deleted □ modified: 
MOTE: 'Domestic priority Information. This information includes, the application number, the filing date, the 
status Ondudlng patent number tf available), and relationship of each application for wNch a benefit 
is chimed under 35 U.S.C. 1 19(6). 120, 121, or 3BSfel Providing this information in the application data 
sheet constitutes the specific reference required by 35 U.S.C. 119(e) or 120, and $ 1.T8fa)(2) or 
§ 1.78(a)(4), and need not otherwise be made part of the specification." 37 C.F.R § 1-76Q>X5). 

0 Domestic priority for this application is claimed as follows; 

Q 35 IXS.C § 119(e): Application No.; 60/a73ro-198 

Filed: PMflBtelT il7i 1999 

Status: Afrfrnfonefl 

Relationship: Prmrf fl-fnnal fnr Nnn-Pirovifiional 

Application No.: , — . 

Filed: 

Status: — — — 

Relationship: 

Application No.: 

Filed: 

Status: . 

Relationship: : 

□ 35 US.C. § 365(c): Application No,: 

Filed: 

Status: , — 

Relationship:. _ — „ — . , 

(Supplemental Application Data Sheet [4-1.2]— page 5 of 6) 
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6, Foreign priority information is being □ added □ deleted □ modified: 

NOTE: "Foreign priority information. This informaticn includes the application number, country, and tiling date 
of each foreign application for which priority is claimed, as weli as any foreign application having a filing 
date before that of the application for which priority is claimed. Providing this information in the 
application data sheet constitutes the ciaim for priority as required by 35 US.C. 119(b) and § 1,55(a). a 
37 C.F.R § 1.76(b)(6)- 

□ Foreign priority Is claimed for this application as follows: 

Country: 

Application No.: 

Filing date: , ,. 

Status: 

Foreign application having a filing date before that of the above application for 
which priority is claimed. 

□ None 

□ Country: „ 

Application No.: , .. , 

Filing date: 

Status: - 

7. Assignee Information Is being □ added □ deleted □ modified: 

NOTE: 37 § U$(b) "Assignee Information. This information includes the name (either person or juristic 
entity) and address of the assignee of the entire right, title, and interest in an application. Providing this 
information in the application data sheet doea net substitute for compliance with any requirement of 
part 3 of this chapter to have an assignment recorded by the Office. M 

NOTE-' 37 C.F.R § 1.215(b): "(b) If applicant wants the patent application publication to include assignee 
information, the applicant must include the assignee information on the application transmittal sheet 
or the application data sheet (§ 1.76). Assignee information may not be included on the patent 
application publication unless this information is provided on the application transmittal sheet or 
application data sheet included with the application on filing. Providing this information on the application 
transmittal sheet or the application data sheet does not substitute for compliance with any requirement 
of part 3 of tNs chapter to have an assignment recorded by the Office." 

The assignee(s) of this application Is/are: 

Name of assignee: 

Address of assignee: 



Extent of interest of assignee in appllcati 
Reg. No. 27712 



Tel. No. { 80? 769-4457 RlrhflrH M. KflTTTLAK 

(type or print name of practitioner) 

IBM CORPORATION, IP Law Dept., 972E 

Customer No. 024241 mnn Rw pt Street 

P.O. Address 

Essex Junction, VT 05452 



(Supplemental Application Data Sheet [4-1.2}— page 8 of fi) 
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distmsutioms rvnCTioKS o* random VA*iA»tit$ 
Thus uKiependeace is not required for the first result. Now 



(6-1S) 



d-8 SYNTHETIC SAMPLING (MONT* CARLO TECHNIQUE) 
The methods discussed so far to And the distribution of 
Y = g(Xu • * - » 

have been approximations or too involved to be practical for Urge problem*. 
In this section a very simple and intuitively satisfying method is presented. 
Its only drawbacks are that it requires a digital computer and general para- 
metric results are not obtained, thus limiting the applicability to synthesis. 

It is assumed that Y = g(X x XJ\s known and that the jomt density 

Ar v is known. Now if a sample value of each random variable were 

<^ r x » *u> * - **, ■ • . *« - * a V ^ ue 5 ^ 

could be computed [say y x = g(* u> *« Then If another set of 

sample values were chosen for the random variables (say *i = *» x « - 

then y t -^>*„,.-.,* te ) could be computed. , . v 

If one had ihe lime one could compute many such sample values 01 J. 
The computer actually supplies the speed that makes many such calculations 
possible. There is just one problem. How does the computer select the 
different values of Aj, ^ ...» ^ii? 

If each of the random variables had a uniform distribution between 0 
and I, numbers for each random variable could be chosen from a table of 
random numbers. Actually, computer routines generate pseudorandom 
numbers which may be used- 

Consider the following case. Let the random variables X u — • » x * 
be independent and each uniformly distributed between zero and one, and 
let T «• X u - - - < X*d be a known function, Then the computer program 
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*-l SVNTHeTlC SAMPL1NC (MONK CARLO T£CHN!Ot/tt> 

to compute an approximation to the distribution of Y consists of the follow- 
ing baste steps. 





Generate 20 




random numbers 




and stove as 




*lw » 






Compute 
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Store Vi 




no — 


enough y,'s ^ 





If a plot is desired it may be convenient to plot a usual bar chart or 
histogram. This method simply calls for breaking the range of Y into say 
30 mutually exclusive cells of the same size and plotting vertically the number 
of samples that fell into that ceil. (See Chapter VEtt for more details of 
histograms.) 

The case of uniformly distributed variables was considered. Now let X % 
have a distribution function F AV To obtain a random sample of X % the follow- 
ing procedure may be used. Select a random sample of U which is uniformly 
distributed between 0 and 1. Call this random sample v % . Then F x ~ x (uDi* 
the random sample of X (see Example 6-7). 

For example* suppose that A' is uniformly distributed between 10 and 20. 
Then 
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As another example suppose that X is normally distributed with mean O 
and variance 1. Then a random value of .250 for U would correspond to 
— .67 for a value of X. This resuH follows from a table of the normal distribu- 
tion function. Practically, most computers automatically generate standard 
normal random variables. If this is the case then random sample of a normal 
random variable Y with mean $i and a variance o- may be generated from a 
standard normal random variable X by recalling that 



or 



o* 



(6-19) 



Equation 6-19 can be easily checked by finding the characteristic function 

The only diflercnoe when the random variables are dependent is that the 
dependence must be taken into account when the random samples are gener- 
ated. We assume that the dependence is expressed In terms of conditional 
distributions. If this is not the case the joint distributions can always be 
reduced to the required conditional distributions. 

For illustration consider three dependent random variables: X lt X u and 
We first generate a random sample of X x by the same methods discussed 
above. Call this sample * a - We next generate a random sample of X % using 
ijc »» , by the same method used before. Call this sample Then we 
use* iX to generate x^. Thus nothing changes except that the 

conditionaTdkstnbuoon functions are used in generating the random sampleS- 
With a fast digital computer thousands of simulations can be run in 
reasonable limes, Monte Carlo solutions often involve 10,000 or more 
simulations. An example is given in the next chapter. 
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6-9 SUMMARY 

The purpose of thi» chapter was to consider the important engineering 
problem of finding the distribution of Y « g(X l% X*, - . . t where the 
distribution of the X 4 *s is known. 

First the problem of Y — was considered and this problem was 
solved by (6-2)* Examples were given to illustrate its application. 

Next r =■= was considered, and it wa$ shown that in the case 

of independent random variables the solution involved convolution of the 
density functions or multiplication of the characteristic functions. The case 
where Y is a linear combination of the Jr/s was shown to be only a slight 
generalization of this problem. 

The general problem was then considered and although a general method 
of solution wa* outlined, the difficulty of solution was illustrated by an 
example, and approximations were suggested. 

Two approximations for the general solution were described. First a 
Taylor series approximation, moments, and the central limit theorem were 
used. Then a Monte Carlo method wag suggested. 

References Bl p D3, D4» and F2 provide additional reading. 



6-10 PROBLEMS 

2* JThas a normal density function with mean 1 and variance 2. 

Find the density of Y. 
2. The power P dissipated in a resistor is P = Assume A -> 2 and / i* a 
random variable with a normal density. 

Find the density function of P. 
a. The output of a fun wave rectifier is Y - \Xl Fmd the density function of Y 

when JThas a uniform density from -1 to *M. 
4. The output of a square law detector Is 

Y^aX 2 , o>0. 
Find the density function of Y in terms of the density function f x ofX. 

159 
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7*1 INTRODUCTION TO TOLSRANCE STUDIES 

is usually a good assumption. When the standard deviation of the pans are 
known, then the mean and variance of the output can be computed using the 
approximation developed in the last chapter. 

'Again something must be assumed to describe the output distribution and 
the probability of being out of tolerance. It is suggested that if enough 
variables are involved and the function is approximately a linear combination, 
chen a normal density may be assumed. Then the probability of being outside 
Tolerance limits can be computed^ 

We now show an example of a tolerance problem which is solved using a 
Monte Carlo approach. 

Example 7-2 

This example* a simplification of a problem that actually occurred in 
practice, was worked by the method shown, and the results actually obtained 
in manufacture corresponded with the theoretical results. 

The simplified version is shown in Figure 7-4. The bar will fit within the 
bracket if Y m < Y % and there will be interference (or no fit) if Y* < Y*. In 
the actual case Y* and Y b involved 41 dimensions and the configuration was 
more complicated than simply a sum of lengths. Actually the equations for 
Y a and Y b involved arcs and angles, thus various trigonometric functions 
were involved- 

The problem was solved by finding the probability distribution of Z = 
Y b — Yc by Monte Carlo sampling. Note that if Z > 0 there is no problem, 
while if Z < 0 there will be no fit and the ports cannot be assembled. 

The problem arose because just as production and assembly were about 
to start it was discovered that interference was possible. (As customary, the 
intent was to tighten the tolerances on each part until no interference is 
possible at the worst case, but the designer made a mistake in his worst case 
calculations*) Then the question was, what is the probahJliiy of inter* 
ference? If it is low enough then it would be better to have a few that would 
not fit, rather than wait and spend the extra money to redesign some of the 
parts. 

To find the probability density of Z via a Monte Carlo technique, one 
must have Z = g(*i> - - • . *ti) and must know the joint distribution of 
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v The function g was found from the drawings using trigonom- 

««V A part of the equation is shown in Figure 7-5 simply to illustrate the 

/■Jc/r* - Vf* - (S - + Itf - >Q ~ *uP ~ *" 



2(K+ »1 



a 




- ain-? I tan-' I ( g _ ^ -7! + y; \ 5; / 

(V? r )y ~ - y ]» f Iff - >o - ffuf - ^f^A 

* \ — ** + W) "~ ' 

Figure 7-5 

form. The various dimensions were assumed to be l^epen<lcm and equally 
likely between their upper and lower tolerance limits. The result or WW 
simulations is shown in Figure 7-6. t > 

Note thai the results appear nearly normal and that interference occurred 
71 times ia 8000 ambulations. Based on these results it was decided to produce 
units without a design change and to rebuild those few on which 
did occur. The results of the actual assembly operation corresponded very 
well with the prediction that 71/8000 would noc fit. 

miliary 

Tolerances of parameters can be combined by Monte Carlo methods or by 
Taylor series aptttWiatiOiV Both produce the probability of bemg outside 
certain limits and are less conservative than deterministic totennce studies. 
Note that the analysis described above is the basis for deciding which 
tolerances to assign to what parts. 

172 
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